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ABSTRACT 


An  atrcroA  configuration  for  a  dock  launched  intercept  mission  was  investigated  as  part  of  an  on¬ 
going  wavender  study  by  the  Naval  Postgraduate  School  and  the  NASA  Ames  Research  Center  The  mission 
requiiements  for  the  earner-launched  and  necoverod  airctaft  iiKludcd  Mach  6  cruise  out  to  a  I  (XX)  nautical  mile 
combat  radius  and  20  minutes  of  combat  followed  by  return  to  the  earner.  A  conical-flow  wavender  served  as 
the  starting  point  for  the  aircraft  configuration  A  hyxlrocaibon  scramjet  was  inu^ratcd  with  the  wavender 
body.  The  aft  end  of  the  wavender  was  faired  to  decrease  the  base  area  thereby  reducing  the  transonic  base 
drag.  A  numcncal  optimization  was  then  completed  to  maximize  the  product  of  LTD  (lift  to  drag  ratio)  and 
(specific  impulse).  V'oriables  fur  liie  upumization  included  the  cone  shock  angle  (used  to  denve  the  conical  flow 
wavender)  arxl  the  geometry'  of  both  the  wavender  body  and  the  integrated  propulsion  system.  The  vehicle  was 
constrained  to  a  minimum  volume  of 3240  cubic  foot,  a  maximum  span  of  60  feet  and  a  fixed  length  of  60  feet. 
The  integrated  propulsion  system  was  constrained  to  produce  a  minimum  contraction  ratio  of  12.0  and 
assurance  that  the  cowl  ^wck  was  within  an  acceptable  distance  of  the  inlet  shoulder  of  the  combustor.  The 
optimum  configuration  met  or  exceeded  all  constraints.  L/D  comparisons  vxtc  made  bc^vtx3^  the  integrated 
configuration  (i.e.,  the  subject  of  this  study),  pure  Mach  6  optimizod  waveriders  and  histoncal  trends. 
Additionally,  model  design,  test  media  and  test  parameter  selection  were  studied  for  planned  low'  speed  wirxl 
and  water  tunnel  tests  as  well  as  performance  predictions  ftir  the  planned  wind  tunnel  tests. 
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I.  INTRODUCTION 


The  quest  for  speed  has  been  a  continuous  theme  in  aeronautical  history  The  reasons 
for  this  quest  have  been  many  and  varied,  not  the  least  of  which  has  been  the  pure  thrill 
and  excitement  of  doing  something  in  the  present  that  could  not  be  done  in  the  past 
Engineers,  manufacturers,  pilots  and  spectators  have  all  shared  in  the  excitement  of 
increasing  aircraft  performance  Of  course,  the  excitement  of  going  faster  does  not,  in 
itself,  justify  the  cost  As  noted  by  Kuchemann,  there  is  a  social  aspect  of  aircraft  design 
to  consider  [Ref  I  p.  3]  Going  faster  must  have  a  value  to  society  which  justifies  the 
cost.  Increases  in  speed  have  historically  been  worthwhile  because  they  increased  the 
convenience  and  decreased  the  time  to  travel  a  given  distance,  thereby  increasing  society’s 
productivity.  An  increase  in  any  figure-of-merit  of  aircraft  performance  must  be  balanced 
by  other  factors,  such  as  productivity,  cost,  safety  and  environmental  issues,  in  order  to 
determine  whether  or  not  the  increase  is  worthwhile. 

Kuchemann  envisioned  the  peaceful  future  of  our  planet  as  a  "global  village"  wherein 
which  people  could  reach  all  of  the  world's  major  population  centers  from  any  location  in 
only  a  few  hours.  Kuchemann  argued  that  this  reduction  in  travel  time  was  necessary  to 
overcome  the  reluctance  of  people  to  travel  great  distances  and  that  bringing  about  this 
"global  village"  by  reducing  travel  times  should  be  one  of  the  social  goals  of  aviation. 
Increasing  the  speed  of  aircraft  may  be  the  only  way  to  achieve  Kuchemann's  vision.  From 
our  present  status  in  aviation,  hypersonic  flight  is  the  next  logical  step  in  the  aircraft 
velocity  spectrum.  One  possibility  for  efficient  hypersonic  flight  is  waverider 
configurations.  [Ref  l:pp.  448-451] 
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To  achieve  his  "global  village",  Kuchemann  stated  that  transport  aircraft  of  the  future 
should  be  capable  of  traveling  a  "global  range"  (half  the  earth's  circumference)  m  two 
hours  To  this  end,  he  defined  a  figure-of-merit  by  which  to  compare  various  aircraft 
types,  as  the  fraction  of  the  "global  range"  an  aircraft  was  capable  ot  traveling  m  two 
hours  Figure  11  shows  a  graphical  representation  ot  Kuchemann's  comparison 
represents  the  fraction  of  the  "global  range"  achievable  in  two  hours  by  the  different  types 
of  aircraft  Three  general  aircraft  types  are  considered,  swept,  slender  and  waveriders 
The  aircraft  type  is  basically  defined  by  its  span  to  length  ratio,  s/f  cruise  Mach  number, 

Mq,  and  body/wing  shape,  swept,  slender  or  waverider  Three  additional  lines  ot  interest 
subdividing  the  figure  are  the  Mq=I  line,  the  s/l-0  2  line  and  the  Ps/(""l  line  Span  to 

length  of  less  than  0  2  leads  to  aircraft  generally  inadequate  on  the  airfield  (Ret  1  p  44Q] 
Ps/(=l  is  where  delta  wings  have  a  "nominally  sonic  leading  edge  "  The  figure  shows  that 

waveriders  are  best  suited  for  the  task  of  bringing  the  time  for  traveling  "global  range" 
down  to  a  few  hours  [Ref  1  pp  448-451] 

Waverider  technology  has  a  wide  range  of  possible  applications,  including,  but  not 
limited  to,  the  well-publicized  hypersonic  transports  of  the  future  The  L/D  performance 
of  waveriders  makes  them  a  logical  starting  point  for  improving  the  performance  ot 
Maneuvering  Reentry  Vehicles  (MaRV)  and  extending  the  range  of  Hypersonic  Glide 
Vehicles  (HGV)  [Ref  2:pp  1-5]  Waverider  characteristics  also  lend  themselves  to 
several  space  applications  such  as  aero-assisted  plane  changes  or  orbit  modifications  and 
atmospheric  long  range  cruising  planetary  probes  [Ref  3.p.  9],  Additionally,  waveriders 
show  promise  as  baseline  configurations  for  single-stage-to-orbit  vehicles.  Their 
potentially  superior  L/D  performance  when  operating  at  design  Mach  number  compared  to 
conventional  wing/body  configurations  makes  waveriders  logical  candidates  for  the  next 
generation  of  hypersonic  aircraft. 
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Figure  1.1  Breguet  Ranges  of  Various  Aircraft  Types.  From  Ref.  1. 


A.  HISTORICAL  PERSPECTIVE 

The  history  of  waveriders  apparently  began  in  1959  when  Nonweiler  introduced  the 
idea  of  a  three  dimensional  body  derived  from  the  flowfield  behind  a  planar  shock  [Ref 
4:pp.  521-528],  This  first  waverider  had  a  delta  planform  and  a  caret  shaped  cross 
section.  When  flying  at  the  design  Mach  number,  the  shock  is  attached  to  the  leading 
edges  preventing  spanwise  flow  and  spillage  from  the  lower  to  upper  surface.  The  high 
pressure  flow  behind  the  shock  is  captured  between  the  shock  and  the  lower  surface  of  the 
body  resulting  in  increased  lift.  Figure  1 .2  shows  Nonweiler's  caret  wing. 

In  the  1960's,  1970's  and  1980's,  the  design  of  waveriders  was  extended  to  include 
those  derived  from  other  known  flowfields  such  as  caused  by  right  circular  and  elliptic 
cones.  Figure  1.3  shows  a  conical  flow  waverider.  In  the  late  1980's,  Bowcutt  and 
Anderson  continued  the  progress  in  waverider  development  with  viscous  optimized 
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waveriders  designed  to  maximize  L/D,  Figure  14  shows  the  Bowcutt/ Anderson  optimum 
result  for  a  Mach  6  waverider  design.  [Ref  6;pp.  15-19] 


B.  DESIGN  MOTIVATION 

The  configuration  studied  herein  seeks  to  provide  one  possibility  in  the  next  logical 
step  in  waverider  evolution;  the  design  and  optimization  of  a  mission-specific  waverider. 
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As  Schindel  [Ref  7]  points  out,  optimum  waveriders  are  not  feasible  vehicles  due  to 
several  practical  considerations.  These  include  volume  requirements,  the  effects  of  base 
drag  and  the  practicality  of  the  proposed  vehicle  shape  (construction  and  internal  packing 
requirements).  Additionally,  an  actual  vehicle  based  on  a  waverider  must  be  designed  as 
an  integrated  engine/airframe  unit.  [Ref  7  p  1] 

In  1988  Vanhoy  [Ref  8]  expanded  the  work  of  Bowcutt  into  the  low-speed  regime 
with  wind  tunnel  tests  of  a  Mach  6  viscous  optimized  waverider  Figure  15  shows  the 
waverider  used  in  Vanhoy's  tests  In  selecting  the  geometry  for  his  tests,  Vanhoy  had  to 
address  such  concerns  as  the  ease  of  construction  of  the  selected  shape  and  whether  the 
base  thickness  was  adequate  for  use  with  a  sting  mount  [Ref  8:p.  8]  Such  constraints 
represent  only  initial  considerations  in  trying  to  create  an  actual  vehicle  from  a  waverider 
The  project  described  herein  considers  mission  and  performance  constraints  on  the  full 
scale  vehicle  as  well  as  the  requirements  of  a  low-speed  wind  tunnel  model. 
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RESULTING  WING  AND  SHOCK 

Figure  1.3  Conical  Flow  Waverider.  From  Ref.  5. 
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II.  MISSION  DEFINITION  AND  ANALYSIS 


The  carrier-launched  intercept  mission  was  specifically  chosen  for  the  development  of 
the  waverider  design  discussed  herein.  Waveriders  seem  well  suited  to  this  type  mission 
due  to  their  inherent  advantage  in  low  drag  at  hypersonic  speeds.  Speed  is  an  important 
factor  in  intercept  missions  because  it  translates  into  distance  away  from  what  the 
interceptor  is  protecting.  Waveriders  may  also  bring  a  better  than  average  maneuverability 
at  hypersonic  speeds  to  the  mission. 

A.  MISSION  SPECIFICS 

Building  on  the  SR-Tl  and  X-15  design,  production  and  operations  databases,  a  Mach 
number  of  six  was  chosen  for  the  cruise  speed  as  a  logical,  practical  and  technically 
feasible  step  up  from  current  capabilities.  Carrier  basing  imposes  additional  design 
constraints,  including  takeoff  and  landing  speeds  and  weights,  length  and  span  and  fuel 
type.  It  was  decided  that  the  more  restrictive  constraints  of  carrier  operation  during  the 
initial  design  phase  would  lead  to  greater  flexibility  in  relaxing  constraints,  if  at  a  later  date 
a  similar  land-based  mission  were  desired.  For  example,  design  constraints  such  as  length 
and  weight  ni^ht  be  relaxed  for  a  land-based  waverider.  Figure  2.1  depicts  the  proposed 
waverider  intercept  mission.  Table  2. 1  contains  specific  mission  parameters  chosen  for 
this  design. 
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TABLE  2.1 

MISSION  SUMMARY 


Figure  2.1  Mission  Profile. 


B.  AIRCRAFT  SIZING 

An  aircraft  sizing  study  was  completed  for  the  selected  mission  parameters  to  be  used 
herein  for  baseline  comparison.  The  sizing  was  based  on  historical  trends  following 
Nicolai  [Ref  9]  and  Raymer  [Ref  10],  It  is  important  to  note  that  the  sizing  studies 
performed  for  this  project  are  based  on  existing  historical  and  current  data  bases.  The  lack 
of  data  on  hypersonic  aircraft  makes  predictions  of  aircraft  size  requirements  and 
performance  only  somewhat  better  than  an  educated  guess.  It  is,  however,  important  in 
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helping  to  validate  conceptual  results  Sizing  began  with  weight,  length  and  span 
limitations  imposed  by  carrier  basing.  The  sizing  calculations  are  contained  in  Appendix 
A.  A  summary  of  preliminary  sizing  assumptions  and  results  for  the  mission  delineated 
above  are  given  in  Table  2.2.  A  comparison  of  preliminary  sizing  results  with  the 
optimized  configuration  is  made  in  Section  IV. 


TABLE  2.2 

PRELIMINARY  SIZING  SUMMARY 


Length 

60  feet 

Span 

60  feet 

Aspect  Ratio 

1  35 

Wing  Area 

2670  square  feet 

Maximum  takeoff  gross  weight 

100,000  pounds 

Maximum  landing  weight 

50,000  pounds 

Empty  weight  fraction 

0.42 

C  ENGINE  AND  FUEL  SELECTION 

For  Mach  6  cruise,  a  supersonic  combustion  ramjet  (scramjet)  was  selected  as  the 
appropriate  propulsion  system.  Hypersonic  aircraft  require  a  more  careful  integration  of 
the  aircraft  engine  and  body  than  conventional  subsonic  and  supersonic  aircraft  due  to 
their  use  of  the  aircraft  body  as  compression  and  expansion  surfaces  for  the  engine.  For 
the  waveriders  studied,  the  engine  and  body  were  designed  from  the  beginning  to  work 
together  as  a  single  unit.  Compression  for  the  propulsion  system  is  "mixed"  in  that  it  is 
accomplished  by  both  internal  and  external  means.  External  compression  is  accomplished 
by  a  series  of  oblique  shocks  which  are  created  by  ramps  on  the  lower  forebody.  Internal 
compression  is  accomplished  by  a  final  oblique  shock  generated  by  the  leading  edge  of  the 
engine  cowl.  For  flight  below  about  Mach  3  the  scramjet  is  impractical  and  the  use  of 
more  conventional  engines  becomes  necessary.  It  is  envisioned  that  two  turbojets  similar 
to  J-58s  could  be  used  for  takeoff,  landing  and  low  speed  flight.  For  a  100,000  pound 
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aircraft  this  turbojet  propulsion  system  would  provide  a  thrust  to  weight  ratio  at  takeoff  of 
approximately  0  6.  A  longtime  requirement  of  two  engines  for  carrier-based  aircraft  has 
been  standard  practice  in  the  U.  S  Navy  and  is  used  in  the  mission  analysis  and  sizing 
studies  for  this  project.  Table  2  3  contains  a  summary  of  engine  sizing  and  performance 
The  two  primary  fuels  considered  for  this  design  were  hydrogen  and  ethylene. 
Ethylene,  C2H4,  was  chosen  as  the  fuel  for  this  design  for  several  reasons.  One  of  the 
primary  concerns  for  fuel  selection  was  shipboard  handling.  Another  major  factor 
considered  was  the  low  density  and  therefore  higher  volume  requirements  of  hydrogen. 
This  creates  problems  for  both  storage  space  onboard  ship  and  increased  aircraft  internal 
volume  requirements  to  fly  the  same  mission  Existing  size  constraints  for  carrier-based 
aircraft  would  limit  the  range  of  a  hydrogen  powered  waverider.  Finally,  ethylene  is  very 
similar  to  JP-7,  the  fuel  used  by  current  J-58  turbojets.  This  allows  the  entire  flight  to  be 
flown  using  a  single  fuel. 


TABLE  2.3 

ENGINE  SIZING/PERFORMANCE  SUMMARY 


0,55-0.8 

Required  TAV  (cruise) 

0.25 

Specific  Impulse  (hydrocarbon  scramjet) 

1300  seconds 

J58  Max  takeoff  thnjst  (sea  level) 

30,000  pounds 
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ill.  VVAVERIDER  GENERATION  AND  OPTIMIZATION 


The  waveriders  analyzed  for  this  project  were  generated  using  the  Waverider  Code 
developed  by  the  Systems  Analysis  Branch  at  the  NASA  Ames  Research  Center  The 
Waverider  Code  and  its  subroutines  are  a  subset  of  the  Hypersonic  Aircraft  Vehicle 
Optimization  Code  (HAVOC),  also  created  by  the  Systems  Analysis  Branch  The 
computer  used  was  a  Silicon  Graphics  Iris  workstation  The  design  goal  of  this  project 
was  to  achieve  a  hydrocarbon-scramjet-powered  waverider  optimized  for  mission 
performance,  while  keeping  computer  processing  time  to  a  minimum  and  checking  the 
results  enroute  to  solving  for  the  optimum  configuration. 

Due  to  the  much  increased  computer  processing  time  of  the  nose-to-tail  flow  solution 
by  the  Waverider  Code  when  using  hydrocarbon  fuels  (vice  hydrogen),  it  was  decided  to 
begin  by  optimizing  pure  waveriders  and  then  hydrogen-powered  waveriders  to  generate  a 
starting  point  for  a  hydrocarbon-powered  waverider  optimization  A  non-optimized 
waverider  had  to  be  generated  which  met  the  optimization  constraints  to  be  used  as  a 
starting  point  for  optimization  due  to  difficulties  experienced  by  the  optimization  routine 
when  starting  outside  the  defined  constraints.  This  was  accomplished  by  generating  over 
one  hundred  waveriders  to  develop  an  understanding  of  how  each  variable  affected  the 
final  waverider  geometry  and  then  finally  moving  toward  a  set  of  design  parameters  that 
produced  a  waverider  which  met  the  desired  constraints. 

To  save  computer  processing  time  and  be  able  to  more  quickly  determine  an  optimized 
waverider's  suitability,  it  was  decided  to  begin  by  optimizing  a  "pure"  waverider  (no 
engine,  no  base  fairing)  for  maximum  L/D.  The  result  of  this  optimization  would  in  turn 
be  used  as  the  initial  geometry  for  a  waverider  with  a  base  fairing,  optimized  for  maximum 
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L/D  The  resulting  optimized,  faired  waverider  was  used  as  the  starting  point  for 
optimizing  a  hydrogen-scramjet-powered  waverider  The  objective  function  for  this 
powered  waverider  was  changed  from  L/D  to  the  product  of  L/D  and  Ijp  Using  the 
product  of  L/D  and  Igp  resulted  in  tradeoffs  during  the  optimization  process  between 
aerodynamic  and  engine  performance.  Finally,  the  result  of  this  optimization  was  used  as 
the  starting  point  for  a  hydrocarbon-scramjet-powered  waverider  optimized  to  maximize 
the  product  of  L/D  and  Ijp,  Figure  3  1  illustrates  the  design  process  used 
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A.  WAVERIDER  GENERATION 


The  Ames  Waverider  Code  begins  its  computations  by  constructing  a  general  cone- 
derived  waverider  from  a  user-defined  flowfield  as  explained  in  Appendix  B  It  departs 
from  the  resulting  "pure"  waverider  with  the  addition  of  an  integrated  forebody  ramp 
system,  scramjet  engine  and  nozzle  At  the  discretion  of  the  designer,  the  lower  surface  is 
faired  at  the  aft  end  to  reduce  the  area  of  the  base  and  thereby  decrease  the  base  drag  at 
low  supersonic  and  transonic  speeds.  The  program  then  computes  a  nose  to  tail  solution 
of  the  conical  flowfield  under  the  vehicle  by  solving  the  Taylor  -  Maccoll  equation  for 
conical  fiow  The  code  creates  a  data  file  which  defines  the  geometry  This  file  is  stored 
for  later  use  to  generate  a  three  dimensional  visualization  of  the  body 

The  "nose  to  tail"  solution  of  the  conical  fiowfield  includes  a  2-D  analysis  of  the  inlet 
and  nozzle  fiowfield  and  a  quasi  1-D  analysis  inside  the  combustor  Both  analyses 
consider  real  gas  effects  Table  3  1  and  3  2  list  highlights  of  each 

TABLE  3.1 

COMBUSTOR  FLOWFIELD  ANALYSIS  SUMMARV _ 

Quasi  1-D  Analysis _ 

Hydrogen  and  Hydrocarbon  Fuels  (or  mixtures) _ 

Multiple  injector  stations  (Mixing  model  -  Combustor  Efficiency) _ 

Augmented  Pre-burning  capability _ 

Heat  balanced 

-  Ramp/Nozzle  heat  transfer  rate  input  vs.  Mach  no. 

-  Combustor  skin  friction  coefficient  input  vs.  Mach  no 


TABLE  3.2 

_ INLET/NOZZLE  FLOWFIELD  ANALYSIS  SUMMARY 

2-D  planar,  inviscid  analysis _  _ 

Real  gas,  equilibrium  or  frozen  flow  in  nozzle  with  mass  fraction  specified 

Ramp/Cowl  geometry  input  as  a  function  of  Mach  no. _ 

Oblique  shocks,  isentropic  waves  and  sliplines  (contact  surfaces)  computed 
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B.  OPTIMIZATION 


A  numericai  optimization  was  achieved  using  VMA  Engineering’s  Design  Optimization 
Tools  (DOT)  DOT  is  an  optimization  routine  which  varies  user-selected  parameters  to 
find  the  local  maximum  or  minimum  value  of  a  selected  design  objective  function  (i  e  .  L/D 
or  L/D’^Igp)  while  satisfying  multiple  user-defined  constraints.  It  uses  finite  difference 
methods  to  calculate  the  gradients  of  the  objective  and  constraint  functions  DOT  solves 
this  nonlinear  optimization  problem  iteratively. 

Specifically,  DOT  uses  the  Sequential  Linear  Programming  Method  to  implement  the 
Modified  Method  of  Feasible  Directions  algorithm  and  find  an  optimum  solution  (in  this 
case  one  which  maximized  L/D’^Ijp  while  still  meeting  all  of  the  constraints)  This  method 
creates  a  Taylor  series  to  approximate  the  objective  and  constraint  functions  and  uses  the 
resulting  approximation  for  optimization,  instead  of  the  original  nonlinear  functions  This 
allows  for  easier  calculation  of  the  objective  and  constraint  functions  during  the 
optimization  Additionally,  the  objective  and  constraint  gradients  can  be  taken  directly 
from  the  linear  Taylor  Series  expansion  DOT  uses  satisfaction  of  the  Kuhn-Tucker 
Conditions  to  determine  when  it  has  reached  a  solution  Three  additional  considerations 
are  checked  to  prevent  unnecessary  iterations  The  routine  will  end  if  a  maximum  number 
of  iterations  (defined  by  the  user)  is  exceeded  or  if  conflicting  constraints  prevent  a 
feasible  solution  within  20  iterations  The  third  condition  is  a  user  defined  tolerance  which 
sets  a  minimum  required  change  in  the  objective  function  after  each  iteration.  If  this 
condition  is  not  met,  asymptotic  convergence  is  indicated  and  the  optimization  is 
terminated.  [Ref  1  l  :pp.  E-3  -  E-36] 

C.  OPTIMIZATION  METHODOLOGY 

A  specific,  iterative  methodology  for  the  optimizat’on  process  was  selected  to  ensure 
consistency  throughout  the  many  optimization  cycles.  Figure  3.2  illustrates  the  selected 
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methodology  Each  optimization  began  by  defining  the  constraints,  variables  and  objective 
function  for  the  optimization  and  visually  checking  their  assignment  in  the  various  code 
and  data  files.  The  affected  files  were  then  renamed  according  to  run  number  and  saved 
This  resulted  in  time  savings  and  increased  ease  when  trying  to  reproduce  optimization 
results.  The  next  step  in  the  selected  optimization  process  was  to  find  a  geometry  which 
met  the  defined  constraints,  to  act  as  an  initial  starting  point  for  the  optimization  routine 
Experience  in  how  each  variable  effected  the  final  geometry  was  used  to  generate  these 
initial  configurations.  After  running  the  optimization  routine,  the  results  were  analyzed 
This  included  checking  to  ensure  that  the  waverider  met  all  constraints,  had  an 
aesthetically  pleasing  geometry  and  delivered  acceptable  performance.  If  the  optimization 
either  failed  or  delivered  unacceptable  results,  the  variable  or  constraint  responsible  had  to 
be  determined  and  corrected  and  the  process  begun  again. 
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Figure  3.2  Optimization  Flowchart. 


D.  FIXED  GEOMETRY  PARAMETERS 


Several  parameters  were  selected  to  remain  fixed  during  the  optimization  process 
These  parameters  were  derived  from  mission  and  performance  requirements  and  included 
the  vehicle  length,  combustor  design,  fairing  start  point  and  equivalence  ratio 

1.  Vehicle  Length 

Vehicle  length  is  constrained  by  aircraft  carrier  operations  A  sixty  foot  length  was 
selected  as  the  maximum  total  vehicle  length  based  on  current  carrier-based  aircraft. 
Table  3.3  lists  lengths  of  various  carrier-based  aircraft  upon  which  the  length  selection  was 
made  [Ref  12]. 

TABLE  3.3 


SIZE  Cl 

3MPARISONOFC/ 

IRRIER'BASED  AIRCRAFT 

AIRCRAFT 

MISSION 

SPAN 

LENGTH 

A-6  Intruder 

all  weather  attack 

53  ft  0  in 

54  ft  9  in 

F-14  Tomcat 

fighter 

64  ft  2  in 

62  ft  8  in 

E-2  Hawkeye 

early 

waming/control 

80  ft  7  in 

57  ft  7  in 

S-3  Viking 

ASW 

68  ft  8  in 

53  ft  4  in 

F-4  Phantom  II 

38  ft  1  in 

63  ft  0  in 

F-18  Hornet 

37  ft  6  in 

56  ft  0  in 

2.  Combustor 

The  SCRAMJET  combustor  (powerplant)  uses  a  three  inch  throat  height  and  five 
foot  combustor  length.  Combustor  width  was  allowed  to  vary  during  the  optimization,  as 
discussed  later.  Maximum  combustor  efficiency  was  assumed  to  be  0.95  for  this  design. 
Combustor  efficiency  is  the  fraction  of  injected  fuel  to  burned  fuel.  Tables  3.4  and  3  5  list 
locations  of  expansions  and  injectors  inside  the  combustor.  Figure  3.3  shows  combustor 
geometry.  Internal  locations  in  the  combustor  are  given  by  station  numbers  which  are 
multiples  of  throat  height  (i.e.,  station  number  1  is  located  3  inches  from  the  throat). 
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TABLE  3.4 

COMBUSTOR  EXPANSION 


TABLE  3.5 

INJECTOR  GEOMETRY 


Figure  3.3  Fixed  Combustor  Geometry. 
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3.  Equivalence  Ratio 

The  equivalence  ratio  is  a  measurement  of  how  close  to  stoichiometric  a  reaction  is, 

as  defined  by  the  following  equation 

^  _  fuel  I  oxidizer 
^  ~  (  fuel  I  oxidizer  )s,auM.om^c 

(|>  =  1  indicates  a  stoichiometric  reaction,  (j)  <  1  indicates  a  fuel  lean  (air  rich) 
mixture  and  4)  >  1  indicates  a  fuel  rich  (air  lean)  mixture.  A  stoichiometric  reaction  is  one 
where  all  of  the  oxygen  atoms  in  the  oxidizer  react  chemically  to  appear  in  the  products 
[Ref  13:pp.  63-69]  A  simplified  example  of  the  reaction  seen  with  the  combustion  of 
ethylene  is  given  by: 

CjH^  +  XOj+X(3J6)Nj  ->aCO^+bH^O  +  {3J6X)Nj 
Reactants  Products 

For  a  balanced  equation,  X,  a  and  b  must  be  defined  such  that: 

Cj//,  +  30j  +  3(3.76)yVj  -►  2COj  +  2HiO  + 1  l.ZSAj 

The  air  to  fuel  (AFR)  ratio  necessary  to  provide  stoichiometric  combustion  is 
determined  by  the  following  equation: 

AFR  =  =  (^'*'3.76(3))  _  14  28molcs  air/  mole  fuel 

On  a  weight  basis: 

AFR=  14.28(28.95)/ 28=  I4.76units  air/ unit  fuel 

An  equivalence  ratio  of  0.6  was  used  for  the  optimization  of  the  hydrocarbon- 
powered  waverider.  This  value  was  chosen,  after  experimenting  with  different  values,  as 
the  one  which  prevented  thermal  choking  in  the  combustor  and  also  resulted  in  the 
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optimizer  driving  the  combustor  width  to  an  acceptable  value.  Thermal  choking  is  defined 
as  the  supersonic  flow  in  the  combustor  being  reduced  to  Mach  1  by  heat  addition  [Ref 
14:pp.  77-85] 

4.  Fairing  Start  Point 

The  parabolic  fairing  designed  by  Pessin  [Ref  15]  was  added  to  the  lower  aft  end 
of  the  waverider  (see  Figure  3  .4).  The  fairing  was  used  to  reduce  the  base  area,  thereby 
reducing  the  base  drag  for  increased  transonic  and  low  supersonic  performance  In  doing 
this,  however,  several  penalties  were  suffered.  Among  these  were  a  reduction  in  internal 
volume  and  a  loss  of  L/D  performance.  There  were  two  basic  considerations  in  selecting 
the  fairing  start  point.  A  fairing  start  located  farther  forward  from  the  base  resulted  in  a 
sma'ler  initial  angle  for  the  fairing,  decreasing  the  effects  of  flow  separation.  The  forward 
point  was  limited,  however,  by  the  requirement  of  the  fairing  to  be  concave  upward.  If  the 
fairing  began  too  far  forward  the  result  was  a  convex  fairing,  resulting  in  increased  drag. 
Additionally,  the  cruise  L/D  performance  went  down  as  the  length  of  the  fairing  was 
increased.  This  made  a  tradeoff  between  initial  fairing  angle  and  fairing  length  necessary 
Since  the  optimization  routine  does  not  consider  the  increased  problem  caused  by  base 
area  at  transonic  and  low  supersonic  speeds,  it  was  not  possible  to  use  the  fairing  start 
point  as  an  optimization  variable.  Thus,  experimentation  with  different  values  using  the 
Waverider  Code  was  required  to  determine  an  acceptable  fairing  start  point.  This 
experimentation  led  to  a  choice  of  start  location  at  85%  body  length  or  5 1  feet  from  the 
aircraft  nose. 

5.  Assumptions 

Overall  vehicle  density  was  assumed  to  be  25  Ibs/ft^.  A  standard  atmosphere  was 
used  for  ail  calculations  and  the  ratio  of  specific  heats  was  assumed  to  be  1.4. 


22 


E.  OPTIMIZATION  CONSTRAINTS 

Constraints  were  placed  on  the  optimization  of  the  waverider  in  order  to  meet 
projected  mission  and  engine  requirements  and  to  produce  an  aesthetically  pleasing  body 
shape.  Specific  constraints  used  during  the  optimization  included:  internal  body  volume, 
span  to  length  ratio,  contraction  ratio,  maximum  leading  edge  temperature,  level 
unaccelerated  cruise,  and  a  requirement  that  the  cowl  shock  hit  the  combustor  inlet 
shoulder 

1.  Volume 

Required  volume  was  based  on  the  preliminary  mission  analysis  and  sizing  studies 
as  discussed  in  Section  II.  The  actual  constraint  used  was  a  volume  factor  of  15%  of  the 
body  length  cubed.  For  a  60  foot  long  waverider  this  gives  a  minimum  volume  of  3240 
cubic  feet.  This  results  in  an  81,000  pound  aircraft  at  the  assumed  density  of  25  Ibs/ft^. 

2.  Span  to  Length  Ratio 

A  maximum  S/L  of  1.0  was  chosen  primarily  to  give  a  balanced  appearance. 
Additionally,  this  leads  to  a  reasonable  maximum  wingspan  at  the  trailing  edge  of  60  feet 
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for  carrier  based  operations.  Table  3.1  lists  spans  of  various  carrier-based  aircraft  for 
comparison. 

3.  Contraction  Ratio 

Contraction  Ratio  is  a  measure  of  external  compression  based  on  the  size  of  the 

captured  streamtube  and  throat  height: 

_  captured  freestream  area 
inlet  area 

Contraction  Ratio  was  constrained  to  a  minimum  of  12.0  for  this  optimization. 
Figure  3.5  shows  the  external  compression  shocks  and  contraction  ratio  geometry. 


Figure  3.5  Contraction  Ratio  Geometry. 


4.  Unaccelerated  Level  Flight  Requirements 

To  achieve  unaccelerated,  level  flight  during  the  cruise  phase  of  the  mission,  the 
forces  acting  on  the  aircraft  must  be  balanced.  Normally  these  forces  are  lift,  drag,  thrust 


and  weight.  For  unaccelerated  flight,  thrust  (T)  must  equal  drag  (D)  Similarly,  for  level 
flight,  lift  (L)  must  equal  weight  (W)  At  the  speed  and  altitude  considered  for  this 
mission,  however,  another  force  becomes  significant  It  is  the  centrifugal  force  generated 
by  flying  a  flight  path  which  accounts  for  the  curvature  of  the  earth  This  centrifugal  force 
acts  in  concert  with  the  lift  force  to  balance  the  aircraft's  weight,  thereby  reducing  the 
amount  of  lift  required  for  level  flight  Weight  At  Altitude  (Wgit)  is  defined  as  the  aircraft 
weight  minus  the  centrifugal  force  Lift  was  required  to  equal  during  cruise  These 
requirements  are  summarized  below  and  depicted  in  Figure  3  6. 

T  =  D 


where;  T=thrust 
D=drag 
W=weight 
F=centrifijgal  force 
L=lift 

Wait=weight  at  altitude 

S.  Leading  Edge  Temperature 

Maximum  leading  edge  temperature  was  constrained  to  3250  °R. 
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6.  Shock  on  Shoulder 

During  an  early  optimization  run  of  a  hydrogen  powered  waverider,  the 
optimization  process  was  requiring  the  cowl  to  be  too  far  forward  in  order  to  increase  the 
objective  function.  This  resulted  in  increased  compression,  thereby  increasing  engine 
performance,  by  using  reflected  shocks  inside  the  cowl  prior  to  the  throat.  On  an  actual 
aircraft  boundary  layer  effects  behind  these  reflected  shocks  would  probably  lead  to  engine 
unstart.  In  order  to  prevent  this  process,  a  constraint  was  introduced  into  the  optimization 
routine  which  required  the  cowl  shock  to  fall  on  the  inlet  shoulder.  This  constraint  was 
met  by  calculating  a  "Shock  on  Ramp"  ratio.  and  Xg  are  the  locations  in  percent  body 
length  of  the  cowl  shock  and  ramp  intercept  and  inlet  shoulder,  respectively.  A  tolerance 
of  approximately  1  inch  prior  to  the  inlet  shoulder  was  allowed  in  the  optimization.  For  a 
body  length  of  60  feet,  this  leads  to  a  Shock  on  Ramp  ratio  of  0.998.  Figure  3.7  shows 
this  "shock  on  ramp"  geometry. 
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Figure  3.7  Shock  on  Ramp  Geometry. 


F.  OPTIMIZATION  VARIABLES 

Optimization  variables  included  parameters  which  defined  the  waverider  body  and 
propulsion  system  geometry  and  the  flowfield  from  which  the  waverider  was  derived. 

1.  Generating  Curve 

Input  geometry  variables  include  the  first  6  coefficients  which  define  the  freestream 
trailing  edge  or  generating  curve  (shown  below).  C7  and  Cg  were  set  to  zero  for  this 
design.  The  generating  curve  is  described  by  Figure  3.8  Table  3  6  lists  coefficient  limits 
for  the  optimization. 

.V  =  C,lan^,  +C,y  +  10Cjy*  +100C^K^  +IOOOC5r*  +IO,OOOCtK‘  +CTCOs(1.570796C,r/tan6».) 
where:  65  =  shock  angle  of  the  conical  flowfield 
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Figure  3.8  Generating  Curve  Geometry  (Aft  View). 


TABLE  3.6 


GENERATING  CUF 

IVE  COEFFICIENTS 

lower  limit 

initial  value 

c, 

0.1 

0.98 

0  779 

C2 

0.0 

0.95 

0.047 

C3 

0.0 

1.0 

0027 

C4 

-0.5 

2.0 

0.044 

_ 

-0.5 

2.0 

0.054 

_ 

0.0 

2.0 

0.014 

2.  Propulsion  System  Geometry 

A  series  of  ramps  on  the  underside  of  the  aircraft  is  used  to  generate  shocks  which 
achieve  the  external  compression  of  the  air  used  by  the  scramjet.  The  geometry  of  these 
"compression  ramps"  is  described  by  the  angle  and  start  locations  of  the  ramps.  Engine 
cowl  geometry  variables  include  the  cowl  start  point  and  length.  The  initial  nozzle  angle 
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was  also  allowed  to  vary  for  the  optimization  Table  3  7  lists  the  limits  and  initial  values 
for  these  variables  used  during  the  optimization  of  the  hydrocarbon  configuration  Figure 
3  9  depicts  propulsion  geometry 


TABLE  3.7 

PROPULSION  SYSTEM  GEOMETRY  VARIABLES 


r— 

lower  limit 

initial  value 

027 

04 

0.33 

2.0 

5.0 

49 

0.42 

0  55 

0.53 

2.0 

70 

6  8 

1  shoulder  location 

0.602 

069 

0  64 

0.55 

065 

061 

0.25 

0.33 

0  3 

1 

20.0 

400 

25.0 

angles  in  degrees 

locations  in  fraction  of  body  length 
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3.  Leading  Edge  Radius 

The  leading  edge  radius  of  the  waverider  was  allowed  to  vary  between  0625  and 
1.0  inches  The  initial  value  for  optimization  was  set  at  0  5  inches  A  0  5  inch  leading 
edge  radius  was  selected  from  previous  waverider  configurations  as  one  which  would 
result  in  a  leading  edge  temperature  of  less  than  3250  °R 

4.  Generating  Flowfield  Shock  Angle 

The  shock  angle  of  the  conical  flowfield  from  which  the  waverider  was  derived  was 
also  allowed  to  vary  during  the  optimization.  The  best  L/D  performance  of  several  of 
Bowcutt's  viscous  optimized  Mach  6  conical-flow  derived  waveriders  occurred  at  a  shock 
angle  of  12  degrees  [Ref  16;p  1 1]  From  this  information  it  was  decided  to  let  the  shock 
angle  vary  near  this  value.  From  experience  with  prior  waveriders  it  was  noted  that  shock 
angles  above  13  degrees  led  to  poor  L/D  performance  when  the  span  to  length  was 
constrained  to  10  and  shock  angles  below  approximately  11  degrees  resulted  in 
unacceptably  low  internal  volume.  Hence,  the  upper  and  lower  limits  were  set  at  1 3  and 
11  degrees,  respectively.  The  initial  value  for  the  optimization  calculations  was  12  774 
degrees. 

5.  Fairing  Start  Angie 

The  initial  angle  for  the  parabolic  fairine  was  allowed  to  vary  between  10  and  45 
degrees.  The  initial  angle  for  optimization  was  set  at  20  degrees. 

G.  OBJECTIVE  FUNCTION 

As  discussed  previously,  the  product  of  L/D  and  specific  impulse  was  used  as  the 
objective  function  for  optimization.  The  logic  behind  maximizing  this  value  comes  from 
the  Breguet  Range  Equation  (for  jets)  shown  below  Since  the  Breguet  Range  accounts 
for  propulsive  efficiency,  heat  content  of  the  fuel  and  aerodynamic  performance  (L/D),  it 
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can  be  considered  a  "figure-of-merit  for  the  whole  aircraft"  [Ref.  l  ;pp.  6-7]  it  can  be 
seen  from  the  equation  that  increases  in  either  L/D  or  Igp  lead  to  proportional  increases  in 
range  for  a  given  velocity  and  fuel. 


where; 


L 

R  =  Hn—\n-^ 

D  W, 


/ 


V 


y 


R=/.„y  —  \n 

sp 


tv. 


R  =  range 

H  =  calorific  value  of  the  fuel 
Tip  =  propulsive  efficiency 
WfAVi  =  segment  weight  fraction 
V  =  velocity 
Isp  =  specific  impulse 
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IV.  OPTIMIZATION  RESULTS 


As  discussed  in  Section  III,  the  objective  of  this  study  was  to  optimize  a  hydrocarbon- 
powered  waverider  for  the  product  of  L/D  and  Ijp,  while  meeting  several  size  and 
performance  constraints.  A  solution  for  an  optimum  hydrocarbon-powered  waverider  was 
found  by  the  optimization  routine  which  met  all  of  the  constraint  functions.  At  the 
optimum  solution,  two  of  the  seven  constraint  functions.  Shock  on  Ramp  and  Contraction 
Ratio,  were  active.  An  active  constraint  is  one  in  which  the  optimization  routine  is  forced 
to  "march"  along  the  constraint  function  in  order  to  increase  the  objective  function  and 
satisfy  the  constraint,  rather  than  moving  in  the  direction  of  the  gradient.  Fewer  active 
constraints  increases  the  likelihood  of  solving  for  the  true  maximum  of  the  objective 
function  in  the  local  space  being  considered.  The  resulting  flow  properties,  body 
geometry  and  propulsion  system  geometry  of  this  optimum  solution  for  a  hydrocarbon- 
powered  configuration  are  discussed  below. 


A.  FLOW  PROPERTIES 

Table  4. 1  contains  free  stream  flow  properties  and  flow  properties  at  the  cowl  lip  and 
combustor  inlet  for  the  optimum  configuration's  calculated  cruise  altitude  of  85,226  ft. 


TABLE  4.1 

OPTIMUM  CONFIGURATION  FLOWFIELD  SUMMARY 


MACH 

■JtUjcjM 

TEMP 

(R) 

■ESH 

RHO 

(lb/ft3) 

Free  Stream 

6.00 

0.3188 

400.520 

5886.2 

0.000067 

Cowl  Lip 

4.156 

2.388 

739.238 

5528.753 

1.841 

Combustor 

Inlet 

2.493 

17.9 

1474.563 

4612.627 

0.001018 
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B.  OPTIMUM  BODY  GEOMETRY 


Specific  results  from  the  optimized  hydrocarbon-powered  waverider  body  geometry 
are  discussed  individually  below  and  are  summarized  in  Table  4  3,  As  discussed  in  Section 
III,  shock  angle  and  generating  curve  coefficients  were  varied  during  me  optimization 
process.  The  values  of  these  variables  which  led  to  the  optimum  configuration  and  the 
resulting  body  geometry  are  summarized  below.  Figure  4  1  depicts  the  resulting  optimum 
configuration. 

1.  Generating  Shock  Angle 

The  angle  of  the  conic  shock  which  defines  the  flowfield  from  which  the  resulting 
optimum  hydrocarbon-powered  waverider  is  derived  was  selected  as  12.9892  degrees  by 
the  optimization  routine. 

2.  Generating  Curve 

The  coefficients  of  the  generating  curve  which  defines  the  upper  surface  of  the 
optimized  hydrocarbon-powered  waverider  chosen  by  the  optimization  routine  are  listed  in 
Table  4  2.  It  is  noteworthy  that  all  of  the  coefficients  are  positive  This  leads  to  a  simple 
curvature  of  the  upper  surface  and  would  likely  result  in  less  expensive  manufacturing 
costs  for  an  actual  aircraft. 


TABLE  4.2 


Cl 

C7 

Cii 

C4 

Cs 

C6 

0.760425 

0.049392 

0.029621 

0.047527 

0.57519 

0.014328 

3.  Span 

The  S/L  of  the  optimized  configuration  is  0.9419.  This  was  a  direct  result  of  the 
tradeoff  between  the  contraction  ratio  constraint  and  trying  to  maximize  L/D*Isp.  Trends 
from  the  initial  non-optimized  waveriders  indicated  that,  in  general,  L/D  increased  with 
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increased  S/L.  In  order  to  achieve  this  increase  in  S/L,  however,  the  waverider  is 
generated  lower  (farther  from  the  cone  axis)  in  the  conical  flowfield  This  leads  to  smaller 
contraction  ratios  and  reduced  engine  performance  The  resulting  planform  has  a  span  of 
56.51  ft  and  a  plan  area  of  2266  33  ft^.  From  this.  Aspect  Ratio  is  calculated  below  to  be 
1.41  [Ref  17:p.  158]. 


AJi  = 


(56.51)^ 

2266.33 


1.41 


where;  b  =  wing  span 
S  =  plan  area 

4.  Internal  Volume 

Internal  volume  is  important  because  it  determines  the  amount  of  fuel  capable  of 
being  carried  and,  in  turn,  the  aircraft’s  range.  Volume  Factor  for  the  optimized 
configuration  is  2.075  percent  of  the  body  length  cubed.  This  leads  to  a  volume  of 
4481.04  ft3  and  a  weight  of  1 12,025.93  lbs. 

5.  Leading  Edge  Radius 

The  leading  edge  radius  of  the  optimized  configuration  is  0  4995  inches.  This  value 
is  a  result  of  the  tradeoff  during  the  optimization  between  leading  edge  temperature  and 
drag.  As  the  leading  edge  radius  goes  up,  drag  increases  and  temperature  decreases. 
Similarly,  as  the  leading  edge  gets  sharper,  the  temperature  increases  and  the  drag 
decreases.  It  is  desirable  for  a  waverider  to  have  as  sharp  an  edge  as  possible  in  order  to 
keep  the  shock  attached  to  the  leading  edges. 
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Figure  4.1  Resulting  Optimum  Conflguration. 
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TABLE  4.3 

OPTIMIZED  CONFIGURATION  BODY  GEOMETRY 


Length 

Span 

Plan  Area 

AR 

Volume 

Weight 

Lead  Edge 
Radius 

60  ft 

56  51  ft 

1  41 

1 12.025  9  lbs 

0  4995  in 

C.  OPTIMUM  PROPULSION  SYSTEM  GEOMETRY 

The  engine  geometry  parameters  of  the  optimized  hydrocarbon-powered  configuration 
are  summarized  in  Table  4.4.  Combustor  width  of  the  optimum  configuration  is  10  238  ft 
The  forebody  ramps,  cowl  and  nozzle  all  share  this  width.  The  optimum  configuration's 
propulsion  system  geometry  is  depicted  in  Figure  4  2. 

TABLE  4.4 

OPTIMIZED  CONFIGURATION  PROPULSION  SYSTEM  GEOMETRY 
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D.  OPTIMUM  CONFIGURATION  PERFORMANCE 


As  discussed  in  the  previous  chapter,  several  constraints  and  an  objective  function 
were  placed  on  the  optimization  process.  All  constraints  were  met  by  the  resulting 
optimum  hydrocarbon-powered  configuration.  Contraction  Ratio  and  Shock  on  Shoulder 
Ratio  were  active  at  the  optimum  solution.  Table  4  5  summarizes  the  optimized 
configuration's  cruise  performance  as  computed  by  the  Waverider  Code. 

TABLE  4.5 


OPTIMIZED  CON 

FIGURATION  PERFORMANCE  SUMMARY 

CONSTRAINT 

LIMIT 

VALUE  FOR  OPTIMUM 

Contraction  Ratio 

>12.0 

11.940 

Lift 

=  Weight  at  Altitude  (Wai^) 

Thrust 

=  Drag 

Shock  on  Shoulder  Ratio 

>0.998,  <  1.0 

0.996 

Volume  Factor 

>0.015 

0.021 

Leading  Edge  Temp. 

<  3250  R 

2589.5  R 

Span  to  Length  Ratio 

<  1.0 

0.9419 

Cruise  L/D 

4.176 

E.  COMPARISON  WITH  PRELIMINARY  SIZING  RESULTS 

Table  4  6  compares  the  optimized  configuration  with  initial  sizing  results.  The 
optimum  configuration  weight  is  for  cruise.  The  preliminary  sizing  weights  are  at  takeoff. 


TABLE  4.6 

COMPARISON  WITH  PRELIMINARY  SIZING 


PRELIMINARY  SIZING 

OPTIMUM  CONFIG. 

WEIGHT 

100,000  ( max  carrier) 

1 20,600  (first  order ) 

112,026 

LENGTH 

60  ft 

60  ft 

SPAN 

60  ft 

56  ft  6  in 

PLAN  AREA 

2670  ft2 

2266  ft2 

ASPECT  RATIO 

1.35 

1.41 

WING  LOADING 

45  lbs/ft2 

50  lbs/ft2 
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V.  PERFORMANCE  PREDICTION  AND  COMPARISONS 


Aerodynamic  performance  predictions  were  made  using  the  VORLAX  and  Waverider 
codes  at  NASA  Ames  Research  Center.  The  VORLAX  code  is  a  vortex  lattice  code  that 
uses  Polhamus'  Leading  Edge  Suction  Analogy.  Briefly,  Polhamus'  theory  accounts  for 
flow  reattachment  inboard  of  the  leading  edge  when  calculating  total  lift  [Ref  8:p.  20] 
This  flow  reattachment  is  caused  by  the  normal  force  imparted  by  the  vortex  flow 
reenergizing  the  upper  surface  boundary  layer. 

Lowspeed  performance  predictions  (specifically  Cl)  are  used  in  Section  VI  as  a  basis 
for  wind  tunnel  model  sizing  and  to  predict  forces  measured  in  the  tunnel.  A  comparison 
is  made  between  predicted  performance  of  the  optimized  hydrocarbon-powered 
configuration  of  this  study,  Vanhoy's  lowspeed  tested  Mach  6  waverider  [Ref  8],  one  of 
Bowcutt  and  Anderson's  Mach  6  viscous  optimized  waveriders  [Ref  16]  and  historical 
trends. 

A.  COEFFICIENT  OF  LIFT 

The  coefficient  of  lift.  Cl,  is  needed  to  predict  forces  seen  in  the  wind  tunnel.  The 
first  method  used  in  predicting  Cl  was  low  aspect  ratio  wing  theory  as  given  by  the 
equation  below.  The  resulting  value  multiplied  by  a  (in  radians)  results  in  values  for  Cl 
versus  a. 

=  -AR  =  -(1.41)  =  2.215 
2  2 

Additionally,  as  noted  by  Pope  [Ref  18;p.  287],  the  Cl  curve  for  aspect  ratios  below 
three  can  be  approximated  by; 


— =  0.008  +  0.01 8i4/?  (per  degree) 
da 


Finally,  a  vortex  lattice  analysis  was  performed  on  the  optimized  hydrocarbon- 
powered  configuration  of  this  study  by  the  Systems  Analysis  Branch,  NASA  Ames 
Research  Center  using  the  VORLAX  code.  One  result  of  the  VORLAX  code  is  a 
prediction  of  Cl  at  differing  angles  of  attack  and  Mach  numbers.  Cl  values  were 
predicted  by  the  VORLAX  run  for  Mach  numbers  of  0  2  and  0.235  in  anticipation  of 
planned  lowspeed  wind  tunnel  tests.  Angle  of  attack  for  these  computations  was  varied 
from  0  to  20  degrees. 

The  resulting  Cl  versus  a  values  for  both  Low  Aspect  Ratio  Theory  and  the 
VORLAX  code  are  listed  in  Table  5. 1 .  These  results  are  also  plotted  in  Figure  5.1.  From 
the  figure  it  can  be  seen  that  the  slopes  of  the  curves  are  almost  identical.  The  values  from 
the  VORLAX  code  are  slightly  lower  than  those  from  Low  Aspect  Ratio  Wing  Theory 
This  is  most  likely  a  result  of  the  slightly  negative  mean  camber  distribution  of  the 
waverider  body  as  computed  by  the  VORLAX  code. 

TABLE  5.1 


Cl  VALUES  FROM  LOW  AF 

1  THEORY  AND  VORLAX 

Low  AR  Theory 

Pope 

VORLAX  1 

AOA 

Mach  0.2 

Mach  0.235 

0 

0 

0.008 

-0.07546 

-0.07582 

2.5 

0.09663 

0.07145 

0.00421 

0.00398 

5.0 

0.1933 

0.1349 

0.08414 

0.08405 

7.5 

0.2899 

0.19835 

0.16411 

0.16417 

10.0 

0.3865 

0.2618 

0.24276 

0.24295 

15.0 

0.5798 

0.3887 

0.39534 

0.39580 

20.0 

0.7731 

0.5156 

0.53918 

0.53987 
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B.  CRUISE  PERFORMANCE 

Increased  cruise  performance  makes  waveriders  generally  better  suited  to  hypersonic 
flight  than  other  configurations.  Of  course,  this  advantage  in  cruise  performance  is  based 
on  the  condition  that  the  waverider  be  operating  "on  design";  the  Mach  number  it  is 
designed  for  and  at  zero  angle  of  attack.  When  on  design,  waveuders  have  shown 
impressive  results.  Figure  5.2  shows  a  L/D  comparison  between  three  hypersonic 
configurations  [Ref  18:p.  28]. 

L/D  at  cruise  for  the  optimum  hydrocarbon  powered  configuration  of  this  study  is 
substantially  lower  than  Bowcutt  and  Anderson's  Mach  6  viscous  optimized  waveriders. 
As  shown  by  Figure  5.3,  the  L/D  of  their  best  optimum  Mach  6  waverider  (05=12°)  is 
over  8.0  [Ref  16:p.  28].  The  lower  L/D  performance  (4.17)  of  the  waverider 
configuration  of  this  study  is  directly  attributed  to  the  addition  of  the  fairing  and 
propulsion  system.  Table  5.2  summarizes  this  comparison.  Cl  for  the  configuration  of 
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this  study  was  computed  at  Mach  six  and  zero  angle  of  attack  using  the  VORLAX  code 
The  volumetric  efficiency  is  given  by; 


448I.O4'"  ^ 

- = - =  0.12 

S  2266.33 


where:  V  =  internal  volume 


Sp  =  plan  area 

TABLE  5.2 

CRUISE  PERFORMANCE  COMPARISON 


1 

Bowcutt's  Optimum 

Powered  Optimum 

0.12 

0  12 

Cl 

0  045 

0.0084 

L/D 

8.1 

4.2 

C.  LOW  SPEED  PERFORMANCE  COMPARISONS 

Comparisons  of  low  speed  results  by  Vanhoy  [Ref  8]  indicated  good  agreement 
between  Polhamus'  Theory,  a  Mach  6  viscous  optimized  waverider,  and  a  delta  wing  of 
similar  length  and  span  as  shown  in  Figure  5.4  [Ref  8:p.  41],  Similar  performance  should 
be  expected  by  the  configuration  generated  by  this  study  For  comparison,  the  Cl  values 
computed  by  the  VORLAX  code  have  been  overlaid  on  Figure  5.4. 
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Figure  5.2  Hypersonic  ConHguration  Performance  Comparison.  From  Ref.  19. 
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Figure  5.3  Results  for  a  Series  of  Optimized  Mach  6  Waveriders  From  the 
University  of  Maryland.  From  Ref.  16. 
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Figure  5.4  Comparisons.  From  Ref.  8. 
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VI.  TEST  PREPARATION 


Low  speed  testing  is  important  in  determining  the  viability  of  any  aircraft.  It  is  critical 
for  earner  based  aircraft  due  to  their  low  takeoff  and  approach  speeds.  Wind  and  water 
tunnel  tests  are  planned  for  Summer  1993. 

A.  WIND  TUNNEL  TESTS 

Low  speed  wind  tunnel  tests  will  be  conducted  using  a  1 5  inch  long  aluminum  model 
in  the  32  inch  by  45  inch  wind  tunnel  at  the  Naval  Postgraduate  School.  Low  speed 
testing  is  intended  to  investigate  the  speed  regime  below  approximately  1 50  knots;  typical 
of  carrier-based  aircraft  approach  speeds. 

1.  Wind  Tunnel  Description 

The  wind  tunnel  planned  for  use  is  a  closed  circuit,  single  return,  horizontal  flow 
wind  tunnel.  Power  is  provided  by  a  100  horsepower  electrical  motor  coupled  to  a  three- 
bladed  variable  pitch  fan  via  a  four-speed  truck  transmission.  The  tunnel  has  a  contraction 
ratio  of  10: 1  and  a  test  section  length  of  48  inches.  The  maximum  test  section  velocity  is 
260  ft/sec  with  a  nominal  freestream  turbulence  level  of  0.2%.  Hinged  glass  windows  on 
both  sides  and  a  window  in  the  ceiling  provide  easy  access  and  viewing.  [Ref  20; pp.  598- 
599]  Figure  6. 1  shows  the  planned  installation  of  the  waverider  in  the  wind  tunnel. 
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Figure  6.1  Planned  Installation  of  Optimized  Hydrocarbon-Powered  Waverider  in 
NPS  32x45  inch  Wind  Tunnel.  After  Ref.  20. 


2.  Test  Parameters 

Initial  low  speed  testing  will  gather  force  and  moment  data  from  -10  to  +20  degrees 
angle  of  attack.  A  sting  mounted,  six  component  internal  strain  gage  balance  will  be  used 
to  measure  the  forces  and  moments  experienced  by  the  model  in  the  wind  tunnel.  These 
measurements  will  yield  lift,  drag  and  side  forces  and  pitching,  rolling  and  yawing 
moments  [Ref  18:p.  126]. 

B.  WATER  TUNNEL  TESTS 

Flow  visualization  tests  will  be  performed  in  the  Naval  Postgraduate  School  Flow 
Visualization  Water  Tunnel  Facility.  This  is  a  horizontal,  continuous  flow  tunnel  with  a  1 5 
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inch  wide,  20  inch  high  and  60  inch  long  test  section  Glass  constmction  allows  easy 
viewing  of  the  model  during  tests.  Maximum  water  velocity  in  the  test  section  is 
approximately  one  foot  per  second  with  a  nominal  RMS  freestream  turbulence  level  of  less 
than  one  percent.  Flow  visualization  is  affected  by  dye  injection  at  points  of  interest  on 
the  model.  Six  pressurized  lines  route  the  dyes  from  the  dye  reservoir  to  *he  model  [Ref 

21:pp  1-2] 
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C.  MODEL  PREPARATION 

Model  preparation  included  material  selection,  balance  selection  and  model  sizing  for 
both  the  wind  and  water  tunnel  models 

1.  Wind  Tunnel  Model  sizing 

Factors  in  choosing  the  model  size  included  balance  availability,  wind  tunnel  test 
section  size,  and  cost.  The  test  section  size  of  32  by  45  inches  limits  the  span  of  the 
model  to  approximately  1 8  inches  Appendix  C  contains  model  sizing  calculations 

2.  Balance  Selection 

A  listing  of  available  sting  mounted,  internal  strain  gage  balances  was  obtained  from 
the  balance  shop  at  NASA  Ames  Research  Center.  Several  balance  diameters  were 
available.  Due  to  the  small  model  size  required  by  the  NPS  wind  tunnel,  the  0.75  inch 
diameter  balance  was  selected.  An  iterative  process  led  to  the  selection  of  a  1 5  inch  model 
length  and  a  balance  with  a  total  normal  force  capacity  of  50  pounds.  At  20  degrees  angle 
of  attack,  the  model  is  predicted  to  generate  a  normal  force  of  approximately  40  pounds 
(see  Appendix  C). 

3.  Water  Tunnel  Model  Sizing 

The  major  consideration  in  sizing  the  model  for  the  flow  visualization  testing  was 
flow  blockage  in  the  water  tunnel  test  section.  A  limit  of  5%  flow  blockage  was  used  in 
sizing  the  water  tunnel  model.  An  8  inch  long  model  remains  below  this  limit  up  to  about 
22  degrees  angle  of  attack  (see  Appendix  C). 

4.  Model  Material  and  Construction 

Options  considered  for  construction  of  the  models  included  the  computerized 
milling  of  aluminum,  polyboard,  or  laminated  mahogany  models  or  the  casting  of  epoxy 
resin  models.  Mahogany  was  ruled  out  for  both  wind  and  water  tunnel  tests  due  to  its 
inability  to  hold  the  sharp  leading  edges  required  by  the  models.  The  radius  of  the  leading 
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edges  on  the  wind  tunnel  model  will  be  approximately  1  mil  Additionally,  it  would  be 
difficult  to  seal  the  mahogany  model  for  water  tunnel  testing.  Polyboard  was  ruled  out 
due  to  its  lack  of  stiffness.  Epoxy  resin  suits  itself  to  making  many  copies  of  the  same 
model  since  it  is  formed  from  a  mold  The  expense  of  using  it  to  make  just  one  model, 
however,  was  prohibitive.  Aluminum  was  chosen  for  both  models  because  of  its  rigidity, 
its  ability  to  hold  an  edge  and  its  relatively  low  cost.  The  specific  alloy  chosen  was  7075 
aluminum,  as  recommended  by  the  Ames  machine  shop.  Table  6  1  contains  material 
properties  of  7075  aluminum  alloy  [Ref  22;p.  245] 

TABLE  6.1 


70 

75  ALUMINl 

IM  PROPERTIES 

Property 

Density 

Ultimate 

Tensile 

Strength 

Tensile 

Yield 

Strength 

Young's 

Modulus 

Shear 

Modulus 

Units 

Value 

2.80 

523 

448 

71 

26.9  1 

i 
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VII.  FOLLOW-ON  RESEARCH 


The  design  of  any  new  aircraft  creates  research  opportunities  in  every  engineering 
field.  This  is  especially  true  if  the  aircraft  being  designed  is  on  the  frontier  of  technology. 
Hypersonic  waveriders  fit  this  description  and  work  is  required  in  every  field  to  determine 
if  they  are  viable  candidates  for  the  high  speed  aircraft  of  the  future.  Material  and 
structural  analysis,  control  system  design,  thermal  control,  and  application  studies  are  just 
a  few  of  the  areas  that  come  to  mind.  For  the  waverider  configuration  of  this  study,  the 
existence  of  a  set  of  models  makes  intensive,  hands  on  aerodynamic  analysis  possible  in 
the  near  term.  This  should  add  dramatically  to  existing  waverider  data  bases. 

A.  PLANNED  NEAR-TERM  TESTING 

The  next  step  in  this  joint,  on-going  study  between  the  Naval  Postgraduate  School  and 
the  NASA  Ames  Research  Center  is  the  planned  lowspeed  wind  tunnel  and  water  tunnel 
tests  for  which  the  15  inch  and  8  inch  models  are  currently  being  constructed.  These  tests 
are  intended  to  gather  force  and  moment  data  and  flow  visualization  data,  respectively. 
The  results  of  both  of  these  tests  will  undoubtedly  lead  to  more  opportunities  for  research. 

The  initial  water  tunnel  tests  will  likely  show  whether  or  not  the  position  of  the  vortex 
flows  around  the  leading  edges  of  the  waverider  are  stable.  If  they  are  not,  testing  to 
determine  if  they  can  be  made  stable  can  be  studied. 

At  1 5  inches  in  length,  the  wind  tunnel  model  is  large  enough  to  make  the  drilling  of 
pressure  taps  at  a  later  date  relatively  inexpensive.  The  material  it  is  constructed  of  also 
makes  testing  in  the  higher  subsonic  and  supersonic  speed  regimes  possible. 
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B.  ADDITIONAL  AREAS  FOR  RESEARCH 

Several  specific  questions  arose  during  this  study  that  apply  directly  to  the  waverider 
configuration  designed  herein.  Each  could  be  complete  projects  in  themselves.  Among 
these,  modifications  to  the  Waverider  Code  to  append  tails  and/or  a  cockpit  to  the 
waverider  is  a  possibility.  Heat  transfer,  emissions,  structural  requirements,  internal 
packing  studies  and  refined  mission  performance  predictions  should  be  studied 
Additionally,  optimization  of  waverider  configurations  for  other  specific  missions 
including  single-stage-to-orbit,  aero-assisted  orbital  transfer  vehicles,  missiles  and  various 
aircraft  missions  can  be  studied.  Research  in  each  of  these  areas,  and  many  not  listed,  will 
serve  to  move  waveriders  closer  to  actual  vehicles  Figure  7  1  depicts  a  conceptual 
aircraft  based  on  the  waverider  configuration  generated  by  this  study. 
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Figure  7.1  Hydrocarbon-Powered  Waverider  Aircraft  Configuration, 


VIII.  CONCLUSION 


Arguably,  the  most  valuable  result  of  this  research  project  has  been  the  successful  use 
of  the  Ames  Waverider  Code  to  generate  a  completely  integrated  engine/body  waverider 
configuration  optimized  for  a  specific  mission.  The  success  of  the  Waverider  Code  in 
producing  a  viable  powered  waverider,  with  so  many  constraints  placed  on  size, 
performance  and  appearance,  attests  to  the  code's  robust  nature  and  the  programming 
talent  of  those  who  wrote  it. 

It  is  concluded  that  an  integrated  engine/body  aircraft  configuration,  capable  of 
meeting  specific  mission  requirements  (a  deck-launched  interceptor  in  this  case),  can  be 
derived  from  a  conical-flow  waverider.  As  expected,  the  aerodynamic  performance  of 
such  a  configuration,  when  compared  with  an  optimum  pure  waverider,  is  substantially 
lowered  by  the  addition  of  a  propulsion  system  and  fairing.  Additional  penalties  will 
certainly  be  suffered  with  the  addition  of  a  cockpit,  fin  and  control  surfaces. 

This  project  has  studied  several  considerations  required  to  design,  model  and  test  an 
actual  powered  aircraft  configuration  based  on  a  waverider.  It  is,  however,  only  a  first 
look  at  these  requirements.  It  is  recommended  that  studies  in  the  optimization  of 
waveriders  for  application  to  other  real-world  missions  be  undertaken.  These  studies 
should  be  of  integrated  engine/body  designs  with  a  view  towards  optimizing  a  specific 
desired  mission  performance.  Only  this  type  of  research  will  answer  whether  or  not 
waverider  configurations  can  be  viable  vehicles. 
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APPENDIX  A  -  PRELIMINARY  AIRCRAFT  SIZING 

A  preliminary  sizing  study  was  done  following  Nicolai  [Ref  9]  and  Raymer  [Ref  10] 
Preliminary  sizing  was  based  on  historical  trends  and  current  data  bases  and  was  used  for 
comparison  with  the  optimized  hydrocarbon  configuration  of  this  study  The  basic  result 
of  preliminary  sizing  is  to  estimate  gross  takeoff  weight  and  fuel  weight  required  for  the 
mission 


A.  WING  GEOMETRY 

Wing  geometry  is  predicted  for  the  given  mission  based  on  historical  trends. 

1.  Aspect  Ratio 

One  would  expect  that  a  delta-type  planform  waverider  would  have  a  low  aspect 
ratio.  Based  on  historical  trends  for  jet  fighters  (non-dogfighters),  the  aspect  ratio  is 
estimated  at  1.35  [Ref  10;p.  53]. 

2.  Wing  Area 

Based  on  the  predicted  aspect  ratio  and  assuming  the  span  is  equal  to  60  feet 
(S/Lmax^*  0)’  fhe  wing  area,  S,  is  equal  to  2670  square  feet. 

AR  1.35 


B.  THRUST-TO- WEIGHT  RATIO 
1.  Takeoff 

A  statistical  estimation  of  TAV  for  a  jet  fighter  (non-dogfighter)  gives  a  typical 
required  installed  TAV  of  0.6  [Ref  10:p.79].  Using  two  J-58  turbojets,  each  with  30,000 
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pounds  of  thnjst  at  sea  level  [Ref  9:p.  14-6],  fills  this  requirement  for  a  100,000  pound 
aircraft. 


2.  Cruise 

During  level,  unaccelerated  cruise,  the  TAV  must  equal  the  inverse  of  L/D  [Ref 
10;p.  81],  Cruise  L/D  was  assumed  to  be  4  0  for  this  calculation 


I 

[HD] 


0.25 


C.  TAKEOFF  WEIGHT  ESTIMATION 

Takeoff  weight  is  estimated  using  an  iterative  process  based  on  an  initial  guess  of 
takeoff  weight,  crew  and  payload  weight,  and  calculated  fuel  and  empty  weight  fractions 
as  given  below  [Ref  10:p.  12].  Fuel  fraction  assumes  a  6%  allowance  for  reserve  and 
trapped  fuel.  [Ref  10:p.  23] 


1.06(1 - 


where:  Wq  =  takeoff  weight 
Wcrew  “  crew  weight 
'’^payload  =  payload  weight 
Wf  =  fuel  weight 
Wg  =  empty  weight 

Wx  =  mission  segment  fiiel  fractions  as  calculated  below 
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1.  Initial  TakeolT  Weight  Estimate 


The  initial  weight  guess  for  this  aircraft  was  chosen  to  be  100,000  pounds  based  on 
the  maximum  takeoff  gross  weight  for  existing  catapults  from  Figure  A  1  [Ref  I0:p  90] 


CATAPULT 


Figure  A.l  Catapult  End  Speeds.  From  Ref.  10. 


2.  Empty  Weight  Fraction 

Empty-weight  fraction  is  estimated  to  be  0.5  for  initial  sizing.  This  is  based  on 
historical  trends  and  also  on  catapult  takeoff  and  arrested  recovery  weight  limitations. 
Maximum  catapult  takeoff  weight  is  approximately  100,000  lbs.  From  Figure  A. 2, 
maximum  arrested  recovery  weight  is  approximately  50,000  lbs  [Ref.  10:p.  92].  These 
maximums  are  used  as  a  conservative  estimate  for  initial  sizing.  The  use  of  advanced 
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composites  in  construction  of  the  aircraft  could  result  in  a  16%  lighter  empty  weight  [Ret 
9;p  5-2]  This  results  in  an  empty  weight  fraction  of 


(S0,000)(l-0.16) 

~  100,000 


ENGAGEMENT 
SPEED  (KTS) 


MASK  7  ARRESTING  GEAR 


MAXIMUM  LANDING  WEIGHT,  lo’lb* 


Figure  A.2  Arresting  Gear  Weight  Limits.  From  Ref.  10. 


3.  Mission  Segment  Weight  Fractions 
a.  Warm-up  and  Takeo/f 

Weight  fraction  for  engine  start,  taxi,  and  takeoff  is  estimated  from  historical 
trends  [Ref  10:p  104] 


WiAVo  =  0  97 
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b.  Climb  and  Accelerate 

For  supersonic  aircraft  the  climb  weight  fraction  with  beginning  Mach  of  0  1  is 
approximated  by: 

W 

— !- =  0  99 1  -  0  007  M  -  0  0 1 M - 

where;  M  =  ending  Mach  number 

The  beginning  Mach  number  for  this  design  will  be  based  on  the  catapult  end 
speed  for  the  aircraft  weight.  From  Figure  A.  1,  this  is  approximately  1 30  kts  or  M  =  197 
For  a  cruise  Mach  number  of  6.0  the  weight  fraction  for  acceleration  trom  M  =  0  1  must 
be  divided  by  the  weight  fraction  for  acceleration  from  M  =  0.1  toM  =  0197  [Ref  10:p 
105] 

0.991  -  0.007(6.0)  -  0.01(36)  ^  ^ 

~  0.991 -0.007(.197)-0.01(. 039) 


c.  Cruise 


The  Brequet  range  equation  is  used  for  cruise: 


This  leads  to  cruise  weight  fractions  of 


RC 


where:  R  =  range 

C  =  specific  fuel  consumption 

V  =  velocity 

L/D  =  lift-to-drag  ratio 

(1)  Range.  Range  for  each  cruise  segment  (out  and  back)  is  1000  nm 

(2)  Specific  Fuel  Consumption.  From  Figure  A.3,  an  initial  estimate  of  Ijp 
equal  to  1300  seconds  is  made  for  a  hydrocarbon-powered  scramjet.  Specific  fuel 
consumption,  the  inverse  of  specific  impulse,  is  equal  to  7.7xl0"^/sec. 

(3)  Velocity.  Velocity  is  calculated  assuming  a  standard  atmosphere  at  80,000 

feet. 

Ma 

where:  M  =  cruise  Mach  number 

a  =  speed  of  sound 

At  80,000  feet  the  speed  of  sound  is  968.08  fl/sec.  For  a  Mach  number  of  6.0 
this  gives  a  velocity  of  5808.48  ft/sec. 

(4)  L/D.  L/D  is  assumed  to  be  4.0  for  initial  sizing  calculations.  This  is  a 
conservative  assumption  based  on  trends  developed  during  early  waverider  configurations. 
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Figure  A.3  Specific  Impulse  Trends. 


(5)  Cruise  Out 


^2 


-(6076000/rH‘y^’‘*0 

(5808.4iy»/KC><4.0) 


(6)  Cruise  Back 


^  =  0.818 
^2 


d.  Combat 


0.818 


Combat  fuel  requirements  are  a  function  of  engine  performance  and  combat 
duration  [Ref  10;p.  106],  For  the  specified  mission,  the  combat  duration  is  20  minutes. 


It  is  assumed  for  these  calculations  that  combat  will  be  flown  at  low  Mach  numbers  using 
the  turbojets  for  propulsion  TAV  for  this  segment  is  estimated  to  be  0  5 

H' 

1-1.12x10  "(.5)(I200)  =  0.897 

where.  C  =  specific  fiiel  consumption  ( 1/sec) 

T/W  =  combat  thrust  to  weight  ratio 
d  =  combat  duration  (sec) 

e.  Loiter 

Typical  requirements  for  aircraft  include  the  capability  of  20  minutes  of  loiter 
before  landing.  This  allows  for  such  eventualities  as  missed  approaches  and  fouled 
runways/decks.  Specific  fuel  consumption  for  this  mission  segment  is  based  on  the  J-58 
turbojets  and  is  equal  to  1.72xl0‘^/sec.  Loiter  for  jet-powered  aircraft  is  flown  at 
maximum  L/D.  Maximum  L/D  is  estimated  to  be  8.0  for  these  calculations. 

Ur  lS£  -(UOOiwMI.TZxlO"*) 

=  «.»  =0.975 

where;  E  =  loiter  time 

f.  Descent  for  Landing 


Descent  for  landing  is  estimated  from  historical  trends  [Ref  10:p.  107]. 


g.  Landing  and  Taxi  Back 


Landing  is  estimated  from  historical  trends  [Ref  10:p.  107] 


^=.992 

w. 


3.  Fuel  Fraction 

The  fuel  fraction  is  calculated  from  the  total  mission  weight  fraction  given  by  the 
product  of  the  mission  segment  weight  fractions  calculated  above 

w 

=  (.97)(.595)(.8I8)(.897)(.818)(.975)(.990)(.992)  =  0.332 

^0 

1.06(1 -0.332)  =  0.708 

^0 


4.  Resulting  First  Order  Estimation  of  Takeoff  Weight 

=  200  +  400  +  (.708)(100,000)+.42(100,000)  =  120,600/A.v 
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APPENDIX  B  -  WAVERIDER  CONSTRUCTION 


Waveriders  are  generated  in  a  reverse  fashion  from  most  aerodynamic  bodies  in  that 
the  flowfield  from  which  the  waverider  is  derived  is  known.  Waveriders  can  be  generated 
from  the  flowhelds  caused  by  shapes  other  than  themselves  (wedges,  cones,  etc.)  The 
waveriders  of  this  study  were  generated  from  the  flowfields  caused  by  axisymmetric 
circular  cones  at  zero  angle  of  attack.  The  explanation  of  waverider  construction 
contained  herein  will  use  this  type  of  waverider  as  an  example.  The  process  used  for  the 
waveriders  of  this  project  is  straightforward  and  fairly  standard. 

STEPl.  SOLVE  THE  FLOWFIELD 

The  first  step  in  creating  a  conical  flow  waverider  is  to  solve  for  the  flowfield 
generated  by  a  right  circular  cone  of  known  semi-vertex  angle.  The  Taylor  Maccoll 
equation  [Ref  16:p.  7]  is  used  to  generate  the  solution  for  this  flowfield: 

(r - 1) / 2[l  - 1'/  -  (</K,  / d0)*][2y,  +  {dV,  ld0)cot d%  lde>‘  -  {dV,  / dO)  +  (</K,  ldO)(d*V,  /«/<?»)]  =  0 

where;  =  nondimensional  component  of  flow  velocity  along  a  conical  ray 
0  =  ray  angle  from  axis 
Y  =  ratio  of  specific  heats 

STEP  2.  GENERATE  FREESTREAM  UPPER  SURFACE 

The  generating  curve  (see  Section  III)  is  drawn  longitudinally  through  the  cone 
flowfield  to  "carve"  out  the  waverider's  freestream  upper  surface.  Figure  B.l  shows  a 
general  cone-derived  waverider,  the  generating  cone  and  the  conic  shock. 
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fr««str«ain  aurfaca 
compraaalon  aurfaca 


Figure  B.l  Construction  of  General  Cone-Derived  Waverider.  From  Ref.  16. 


STEP  3.  LOWER  COMPRESSION  SURFACE 

Conic  streamlines  are  traced  back  from  the  intersection  of  the  shock  and  the  upper 
surface  to  generate  the  lower,  compression  surface  of  the  waverider. 


APPENDIX  C  -  MODEL  SIZING  CALCULATIONS 

A.  WIND  TUNNEL  MODEL 

The  two  major  factors  considered  in  sizing  the  wind  tunnel  model  were  the  test  section 
size  and  normal  force  measuring  capability  of  the  sting  mounted  balance.  The  test  section 
size  of  the  NPS  tunnel  (32  x  45  indies)  limits  the  span  of  models  to  approximately  18 
inches.  The  small  size  of  such  a  model  indicated  the  need  to  use  a  0  75  inch  diameter 
balance  (the  smallest  available  from  Ames).  The  0.75  inch  diameter  balances  were 
available  with  normal  force  capabilities  of  50  and  100  pounds.  A  general  rule  of  thumb  for 
force  measurements  is  to  use  approximately  75  percent  of  the  balance  capacity  at  the 
highest  force  measured.  This  will  generally  result  in  the  best  resolution  of  measured 
forces.  The  goal  of  this  sizing  analysis  was  to  find  the  largest  model  possible  which  would 
generate  a  normal  force  of  75  percent  while  at  the  maximum  angle  of  attack  planned  for 
the  tests. 

The  normal  force  expected  in  the  tunnel  is  a  function  of  the  dynamic  pressure,  q,  the 
coefficient  of  lift.  Cl,  and  the  size  of  the  model.  For  best  resolution,  this  force  should  be 
37.5  or  75  pounds  for  the  50  and  100  pound  balances,  respectively.  The  maximum  angle 
of  attack  for  the  initial  tests  will  be  20  degrees.  The  predicted  Cl  value  at  20  degrees 
angle  of  attack  was  calculated  by  the  VORLAX  code  at  NASA  Ames  as  0.54. 

For  a  standard  atmosphere  at  sea  level,  the  speed  of  sound  is  1116.4  ft/sec  and  the 
density  is  2.3769x10"^  slugs/ft^.  The  planned  test  velocity  is  140  kts  or  236.3  ft/sec. 
From  this,  the  dynamic  pressure  is  calculated  as  follows: 
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q  =  ^/7v^  =(.5)(2.3769x  I0^)(236.3)^  =  66.36lbs /ft^ 

where:  p  =  density  (slugs/ft^) 

V  =  velocity  (fVsec) 

From  the  optimized  hydrocarbon-powered  configuration,  the  plan  area,  s,  is  given  by 

s  =  (.674)(/)*(.9416) 


The  lift  force  can  be  calculated  by: 


L  =  C^qs 


where:  Cl  =  coefficient  of  lift 

q  =  dynamic  pressure 
s  =  plan  area 

The  100  pound  balance  was  ruled  out  due  to  the  inability  of  an  18  inch  model  (the 
largest  possible  in  the  tunnel)  to  generate  75  pounds  of  normal  force  (lift)  at  20  degrees 
angle  of  attack  as  shown  below: 

L  =  (. 54 )(66.36)(  1.428)  =  51.2  lbs 

The  model  length  necessary  to  generate  37.5  pounds  of  force  can  be  found  in  a  reverse 
fashion: 


/  = 


37.5 


(.54)(66.36)(.674)(.9416) 


=  1.28  ft  =  15.4  in 
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B.  WATER  TUNNEL  MODEL 

The  major  factor  in  sizing  the  water  tunnel  model  is  flow  blockage  by  the  model  in  the 
water  tunnel  test  section.  A  general  rule  of  thumb  is  to  limit  this  flow  blockage  to  five 
percent  of  the  test  section  area.  The  test  section  area  of  the  NPS  water  tunnel  facility  is 
300  in^.  This  gives  a  desired  maximum  blockage  of  15  in^.  The  blockage  is  a  function  of 
the  model  length  and  angle  of  attack.  For  this  series  of  tests  the  maximum  angle  of  attack 
of  interest  is  approximately  20  degrees.  The  maximum  length  of  the  model  can  be 
calculated  by  knowing  the  plan  area  which  generates  the  maximum  blockage  at  20  degrees 
angle  of  attack  and  is  calculated  below 

The  maximum  blockage  is  given  by: 

(.05)(300iV)  =  15iV 

From  this,  the  plan  area  of  the  model  follows  as; 


15in^ 

sin(20°) 


43.86in' 


The  relations  between  length,  span  and  plan  area  are  known  from  the  optimized 
hydrocarbon-powered  configuration  and  are; 

plan  area  =  (.674)(length)(span) 
span  =  (.9416)(length) 


From  this,  the  model  length  is  computed  as: 

length  = 


43.857 


(.674)(,9416) 


=  8.3in 
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From  the  results  above,  a  length  of  eight  inches  was  selected.  This  length  will  result  in 
a  maximum  angle  of  attack  of  approximately  22  degrees  at  five  percent  flow  blockage. 
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